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ABSTRACT 

The real-space optical depth distribution along the line of sight to the QSO Q1422+231 
is recovered from two HIRES spectra using a modified version of the inversion method 
proposed by Nusser & Haehnelt (1999). The first two moments of the truncated optical 
depth distribution are used to constrain the density fluctuation amplitude of the inter- 
galactic medium (IGM) assuming that the IGM is photoionized by a metagalactic UV 
background and obeys a temperature-density relation. The fluctuation amplitude and 
vq ■ the power law index a of the relation between gas and neutral hydrogen (Hi) density 

f^*i ' are degenerate. The rms of the IGM density at z « 3.25 estimated from the first spec- 

CTN ■ trum is a = \/exp [(1.8 ± 0.27) 2 /a 2 ] — 1, with 1.56 < a < 2 for plausible reionization 

histories. This corresponds to 0.9 5s a 5s 2.1 with a{a = 1.7) = 1.44 ± 0.3. The values 
f~| . obtained from the second spectrum are higher by s=s 20%. If the IGM density traces the 

Oh' dark matter (DM) as suggested by numerical simulations we have measured the fluctu- 

q . ation amplitude of the DM density at an effective Jeans scale of about a hundred to two 

5h ' hundred (comoving) kpc. For CDM-like power spectra the amplitude of dark matter 

t/3 . fluctuations on these small scales depends on the cosmological density parameter Vt. For 

power spectra normalized to reproduce the space density of present-day clusters and 
with a slope parameter of T = 0.21 consistent with the observed galaxy power spec- 
trum, the inferred $7 can be expressed as: fi = 0.61 (a/1. 7) ' (a;,/0.62) for a flat 
5_^J ■ universe, and SI = 0.91 (a/1.7) ' (a;,/0.62) for a A = universe. x } is the effective 

Jeans scale in (comoving) /i _1 Mpc. Based on a suit of detailed mock spectra the \-a 
error is ~ 25%. The estimates increase with increasing T. For the second spectrum we 
obtain 15% lower values. 

Key words: cosmology: theory, observation , dark matter, large-scale structure of the 
Universe — intergalactic medium — quasars: absorption lines 



1 INTRODUCTION 

The Lyman forest in QSO absorption spectra is now generally believed to be caused by absorption by large-scale Hi density 
fluctuations of moderate amplitude in a warm photo-ionized intergalactic medium (IGM) . This is substantially different from 
the old cloud picture for the Lya forest (see Rauch 1998 for a review). It is mainly sustained by cosmological hydrodynamical 
simulations which successfully reproduce the absorption features in QSO spectra (Cen et al. 1994; Petitjean, Miicket & Kates 
1995; Zhang, Anninos & Norman 1995; Hernquist et al. 1996; Miralda-Escude et al. 1996, Theuns et al. 1998). The picture 
had, however, been suggested before these numerical simulations became available (Bond, Szalay & Silk 1988; Bi, Borner & 
Chu 1992). Important results of the simulations are a tight correlation between the Hi and the dark matter distribution (on 
scales larger then the Jeans length of the IGM) and a simple gas temperature-density relation which depends only on the 
reionization history of the Universe (Hui & Gnedin 1997, Haehnelt & Steinmetz 1998). As a consequence, the Lya forest can 
be an important probe for the distribution of dark matter and its evolution over a wide redshift range (c.f. Weinberg 1999 
for a review) on scales smaller than those accessible by the galaxy distribution or the lensing effects of foreground large scale 
structures on background galaxies (c.f. Bartelmann & Schneider 1999). 

For example, Gnedin & Hui (1998) showed that the column density distribution of Lya absorption lines depends on the 
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amplitude of the density fluctutions in the IGM. Croft et al used the flux power spectrum of QSO absorption spectra to 
constrain the DM power spectrum on scales from 2/i _1 to 12/i _1 Mpc (Croft et al. 1998, Croft et al. 1999, Weinberg et al. 
1999). Nusser and Haehnelt (1999) (hereafter NH99) proposed an iterative procedure to obtain the real space density and 
its probability distribution function. Here we apply a modified version of this technique to two HIRES spectra of the QSO 
Q1422+231. We infer the fluctuation amplitude of the DM density on the effective Jeans scale and constrain the normalization 
constant of the optical depth distribution which depends on the ionization rate Tphot, the baryon density fibar and the Hubble 
constant (.4 oc fi 2 ar /i 3 /r p hot). Section 2 describes our modified technique to extract moments of the probability distribution 
(PDF) of the density and presents tests on mock spectra. In section 3 we present the application to the data. Section 4 
discusses implications for the cosmological density parameter. In section 5 we relate our work to previous work and give an 
outlook. Section 6 contains our conclusions. 



2 THEORETICAL BACKGROUND 

2.1 Basic equations 

The optical depth in redshift space due to resonant Lya scattering is related to the Hi density along the line of sight (LOS) 
in real space by 



r{w)=o — — I n HI (x)H[w- x-v p (x),b(x)]dx, (1) 

where ctq is the effective cross section for resonant line scattering, H(z) is the Hubble constant at redshift z, x is the real 
space coordinate (in kms -1 ), His the Voigt profile normalized such that J H. dx — 1, v p (x) is the LOS peculiar velocity, 
and b(x) is the Doppler parameter due to thermal/turbulent broadening. The absorption features in the Lya forest are 
mainly produced by regions of moderate densities where photoheating is the dominant heating source and shock heating 
is not important. Hydrogen in the IGM is highly ionized (Gunn & Peterson 1965, Scheuer 1965 ) and the photoionization 
equilibrium in the expanding IGM establishes a tight correlation between neutral and total hydrogen density. This relation 
can be approximated by a power law n HI oc n^, where the parameter a depends on the reionization history. The possible 
range for a is is 1.56 ^ a ^ 2 with a value close to 2 just after reionization, and decreasing at later times (Hui & Gnedin 
1997). Numerical simulations have shown that the gas density traces the fluctuations of the DM density on scales larger than 

the Jeans length, so that n HI = n HI I -°M — J , where n HI is the Hi density at mean dark matter density. In this relation 

Pdm ( x ) i s the dark matter density smoothed on the Jeans length scale below which thermal pressure becomes important. The 
Jeans length in the linear regime is given by (comoving), 

27rCs n , \ 
x, v = 1 + g 

J ' lln ^/^Gp K ' 



l:i ljTP«) I TTTTTnTI? ) ( i 1 ,o „.wn 7 J ( — ) M P C ' ( 2 ) 



- i/2 / t y /2 ( i.5 ^ ' -' 

0.125^ l v 1.5xl0 4 K y ) l y l + (2-a)/0.' 

where c s is the sound speed. However, in the non-linear regime gas can collapse to scales smaller than this and the Jeans scale 
becomes a somewhat ambiguous quantity. Here we define an effective non-linear Jeans length, x 3 , as the width of a kernel of 
the form [1 + (kxj2-K) 2 ]~ 2 , such that the rms fluctuation amplitude of p DM (x) is the same as that of the unsmoothed dark 
matter density filtered with this kernel (see section 4 for details). On scales larger than the effective Jeans length, equation 
((D) can be written as 

t(w) = A(z) f ( Pdm(x) ) H[w-x-v p {x),b(x)]dx, (3) 

J -oo V Pom J 

with 

A(z) 



0.61 



H(z) ■-' 

/ 300km s _1 Mpc _1 \ /fibarfr 5 



V H(z) J \ 0.02 



I r phot \- / t \ (i+±y 

Vl0- 12 s-V l v 1.5xl0 4 K y / V 4 J ' 

(4) 



where fibar is the baryonic density in terms of the critical density and Tphot is the photoionization rate per hydrogen atom. 
The Doppler parameter in the last equation depends on n HI as b oc n H ~ Q ' . 

NH99 have presented a direct Lucy-type iterative scheme to recover the optical depth and the corresponding mass and 
velocity fields in the LOS from the normalized flux, F = exp(— r). They showed that the density field can be successfully 
recovered below a threshold value above which the corresponding flux saturates. NH99 further used mock spectra extracted 



from an N-body simulation to show that the moments of the mass density can be reasonably well estimated by fitting a 
log-normal PDF to the low density tail of the PDF of the recovered density. In the next section we present a modified method 
for estimating the moments of the density distribution. 

2.2 Moments of the density PDF from the recovered local optical depth 



NH99 provided an estimate for the quantity 

■ P( x ) 
r{x) = A 



(5) 



which we term the local optical depth. It is related to the observed optical depth r by a convolution with a Voigt profile as 
described in equation (|l|). The recovered f is a good approximation to the true field only in regions with f less than a certain 
value. For large optical depths, the recovered f typically underestimates the true field. We define a truncated local optical 
depth ft as f t = f for f < f c , and ft = otherwise. The corresponding truncated moments of f t can be written in terms of 
the density PDF as, 



< T t > = 



A n 



(i + sy ia p(s)ds, 



(6) 



where S = p/p — 1 is the density contrast and S c — (f c /A) — 1. We further define v — [ln(p/p) — pi]/p 2 , where fj,i and 
/X2 are the average and rms values of In (1 + §). NH99 had also shown that the PDF of the DM density smoothed on the 
scale relevant for the Lya forest (the effective Jeans scale) can be reasonably well approximated by a log-normal distribution 
(Bi & Davidsen 1997, cf. Sheth 1998 and Gaztanaga & Croft 1998 for different forms of the PDF). For a log-normal density 
distribution, P(v) = exp(— u 2 /2)/y/2~K, the truncated moments in (pi) can be written as, 



< ft >= 



A n 
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where u c is the value of v corresponding to <5 C . By expressing v in terms of f in (pj) we find 
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(8) 
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We have here used the relation p,\ = —p, 2 /2 which follows from the condition < 8 >= for the log-normal. The moments of 
the truncated optical depth distribution depend on four parameters A. , H%, ot and f c . The parameter f c is chosen such that 
for f < f c the local optical depth does not suffer from the biases introduced in saturated regions. As apparent from equation 
( |10|) there are two basic degeneracies leaving us with two independent parameters: 

ant/2, 



B = lnA- 



ap 2 . 



In terms of these two parameters, the moments of ft can finally be written as 
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(10) 



The first two moments, < f t > and < f t >, are therefore sufficient to determine the parameters B and C. From these we 
can then infer the rms fluctuation amplitude of the DM density and the normalization constant of the optical depth A as 
discussed in section 3.3. The use of truncated moments of the local optical depth has a number of advantages. It is numerically 
more stable than fitting a parametric form of the PDF to the distribution of the recovered density. It allows to introduce a 
cut-off in the form of an optical depth instead of the density which depends on the a priori unknown value of A. Finally, it 
shows the basic degeneracies between a, \i 2 and A. 

In reality the PDF of the mass density smoothed on an effective Jeans scale is expected to be somewhat more positively 
skewed than a log-normal. The relations derived above can easily be generalized to other forms of the density PDF. We were, 
however, not able to constrain the additional skewness parameter of an Edgeworth expansion. The same was found by NH99. 
We therefore used a log-normal PDF for the analysis of the data. We assess the possible biases that arise in the determination 
of C and B using mock spectra in the next section. 



3 APPLICATION TO THE SPECTRUM OF Q1422+231 

3.1 The data 

We had two spectra of Q1422+231 at z — 3.62 available for our analysis. The spectra were observed with the HIRES 
spectrograph by Rauch et. al. (1997) and Songaila & Cowie (1996) on the Keck I telescope. The resolution of the spectra is 



Nusser & Haehnelt 



0.100 



Q 



0.010 



0.001 




100.00 



Figure 1. The curves on the right are the PDFs of ln(l + <5) for an Edgeworth expansion around a log-normal distribution with [13 = 0.4. 
The curves on the left show (1 + <5)P[ln(l + <5)]. The solid curves are the average of 70 PDFs from which the mock spectra of Q1422 were 
generated. The dotted lines correspond to a log- normal PDF with the same value of p>2- 

fa 8km s _1 and the pixel size is 2.5km s^ 1 . We used the wavelength region from 4742 A(Ly/3 at the redshift of the QSO) to 
5618 A(Lya at the redshift of the QSO) to avoid confusion due to an overlap of intervening Lya and Ly/3 absorption. This 
corresponds to a redshift range from z — 2.9 to z = 3.62 and is equivalent to 240/i -1 Mpc (comoving) in an Einstein-de Sitter 
Universe. The two spectra were continuum-fitted independently by the observational groups. The differences in our results 
should thus give an idea of the error due to continuum-fitting. 



3.2 Error estimates and tests of the procedure using mock spectra 

Two sources of errors are important for our analysis: first the intrinsic error caused by the inversion procedure; second the 
cosmic variance error resulting from the use of only one quasar spectrum with limited redshift range. We quantify these errors 
and possible systematic biases by analyzing mock spectra. The mock spectra were produced using a modified version of the 
method suggested by Bi et al. (1992) [see also Bi & Davidsen (1997)]. The mock spectra have the same redshift range, pixel 
size, resolution and S/N as the observed spectra of Q1422. We generated random realizations of a Gaussian density field with a 
linear CDM-like power spectrum. We then assumed a parametric form for the PDF of the non-linear density and obtained the 
non- linear density field along the LOS from the linear field by a rank-ordered mapping (Weinberg 1992). The corresponding 
non- linear peculiar velocity, v pcc , was computed as described in NH99. The optical depth was calculated using equation (J3J) 
with a = 1.7 and b = 16n°' I 15 kms _1 . The factor A was chosen such that the mean flux decrement in the mock spectrum 
matches the observed value. Finally, instrumental broadening as well as photon and pixel (read out) noise was added. 

Based on numerical simulations, we expect the true PDF to be more positively skewed than a log-normal. Our method 
for estimating B and C assumes, however, that the PDF is log-normal. It is therefore important to test the sensitivity of the 
estimated moments to the shape of the true PDF. We generated mock spectra for a log-normal and for PDFs described by 
an Edgeworth expansion around the log-normal, 



P(u) = (2,y^ exp (-u 2 /2) 1 + M („ 3 _ 3 „) 



3! 



(11) 



where the skewness parameter is defined as 113 =< v / L12 >. N-body simulations with Gaussian initial conditions have shown 
this parameter to be positive, independent of the detailed shape of the initial power spectrum (Colombi 1993, NH99). We will 
present detailed results for two different PDFs with (1x2,1x3 = 1.3,0.0) and (1.0,0.4), but we have also tested the technique 
with (0.9, 0.6). In figure hi we show representations of the PDF with [13 = 0.4. The deviations of the PDF of (1 + 5) from a 
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Figure 2. Curves of fP[ln(f )] for the recovered and true local optical depth for mock spectra generated with [13 = 0.4. The solid and 
dashed curves are for the recovered t in real and in redshift space. The dotted line correspond to the mean of the true PDFs of the 70 
mock spectra. The error bars represent the cosmic variance error (rms values of the deviations of the mock spectra from the mean PDF). 



log-normal distribution with the same fi2 are significant only for 8 ;>2. Since only the low density part of the PDF is relevant 
for the truncated moments, the method is expected to produce reasonable estimates of ^2 even if the true PDF is somewhat 
skewed relative to log-normal. We have generated mock spectra for both PDFs with and without peculiar velocities to test the 
effect of redshift distortions. For each case 70 mock spectra (i.e. 70 random realizations for the same density power spectrum) 
were analyzed. 

We have also modified the inversion procedure by introducing a power law relation between the Doppler parameter and 
the local optical depth, b — bof 1-01 ' 2 . This has significantly speeded up the convergence of the inversion. To determine 60 we 
define the function \ 2 — J2 ixds (^? bS — ^7*°) /°f 1 where a t . is the error and F; obs the observed flux in pixel i. F/ cc is the 
flux which corresponds the recovered local optical depth. We choose bo to be the largest value which gives % 2 /pixel = 1 after 
a certain number of iterations. 

In Fig. H the average fP[ln(f)] of the recovered local optical depth of 70 mock spectra is shown (for the case with ^3 = 0.4) 
in real (solid curve ) and redshift space (dashed curve). The dotted line represents the average true PDF. The error bars are 
la cosmic variance errors calculated from the 70 random realizations. The spikes in the PDFs of the recovered local optical 
depth are due to the systematic underestimate in saturated regions. Figure H demonstrates that we have successfully corrected 
for redshift distortions in regions of low optical depths. This correction shifts the location of the prominent spike to lower f 
because the redshift distortions tend to enhance the real space density contrast. We use the location of the spike to choose 
the cut-off f c for the truncated moments. 

We have calculated the first two moments < f t > and < f t 2 > for three cutoff values, f c = 1.0, 1.5 and 2.0. The parameters 
B and C were determined by solving equations (|L0l) with Broyden's method (e.g. Press et al 1992). Figures (0) and (0) show 
the true vs the recovered values of the parameters B and C for four sets of 70 mock spectra (two input PDFs, with and without 
peculiar velocities). 

The tight correlation between the true and recovered values demonstrates robustness of the method. The estimated 
values of C are unbiased in all four cases while B is systematically underestimated (except in the log-normal case with peculiar 
velocities set to zero). We found similar results for the mock spectra produced with(//2 = 0.9 , ^3 = 0.6) 

In tables (|l|) and (0) we list the mean values of (C, B) and the estimated errors determined from the four sets of 70 mock 
spectra. The intrinsic random error introduced in the inversion method is the rms value of the differences between estimated 
and true values from each set of the 70 realisations. The cosmic variance errors are estimated as the rms values of the true 
parameters. For both the skewed and the log-normal PDFs the estimated and true values of C agree well within the la errors. 
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Figure 3. The estimated vs true C and for 70 different random realizations of mock spectra generated with v pcc = 0. 
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Figure 4. The same as the previous figure but for spectra generated with non-vanishing peculiar velocities. 



Table 1. True and estimated parameters for the mock ql422+231 spectra generated with a log-normal PDF of the mass density. Listed 
are estimated parameters from mock spectra with and without peculiar velocities (v pcc = 0). The first entry shows the parameters 
estimated after correcting for redshift distortions, the second before the correction and the third the parameters from spectra generated 
with v poc = 0. The first and second rows are (C, B) and the la error due to the the inversion. The fourth entry lists the true values of the 
parameters and the lcr error due to cosmic variance. The transformation from C and B to A is given by In .4 = B + C 2 /2/a and requires 
a value for a. 





fc = 1.0 
(C,B) 


f c = 1.5 

(C,B) 


f c = 2.0 

(C,B) 


Est. (real space) 
inversion error 


(2.27,-2.04) 
(0.1,0.0.22) 


(2.18,-2.22) 
(0.10,0.14) 


(2.16,-2.23) 
(0.11,0.11) 


Est. (redshift space) 
inversion error 


(2.25,-2.25) 
(0.12,0.15) 


(2.23, -2.68) 
(0.11,0.12) 


(2.23,-2.68) 
(0.12,0.11) 


Est. (v pcc = 0) 

inversion error 


(2.23,-1.70) 
(0.08,0.25) 


(2.20,-1.81) 
(0.08,0.18) 


(2.19,-1.85) 
(0.10,0.15) 


True 

cosmic variance 


(2.27,-1.89) 
(0.12,0.16) 







Table 2. The same as the previous table but for mock spectra generated using a skewed PDF with fig = 0.4. Notice that there are 
two entries for true values in this table. This is because different values of A are needed to match the observed flux decrement for mock 
spectra generated with and without peculiar velocities if (13 7^ . 





f c = 1.0 
(C,B) 


T C = 1.5 
(C,B) 


f c = 2.0 

(C,B) 


Est. (real space) 
inversion error 


(2.07,-1.83) 
(0.09,0.18) 


(2.02,-1.95) 
(0.10,0.15) 


(2.00,-1.97) 
(0.10,0.13) 


Est. (redshift space) 
inversion error 


(2.05,-2.26) 
(0.10,0.16) 


(2.03,-2.29) 
(0.11,0.13) 


(2.03,-2.29) 
(0.10,0.10) 


Est. (v pcc = 0) 

inversionerror 


(2.09,-1.95) 
(0.08,0.19) 


(2.07, -2.00) 
(0.08,0.15) 


(2.06,-2.02) 
(0.09,0.13) 


True (v pec ^ 0) 

cosmic variance 


(2.08,-1.30) 
(0.13,0.19) 






True (upec = 0) 

cosmic variance 


(2.08,-1.67) 
(0.13,0.19) 







Redshift distortions only moderately affect the estimated C despite substantial differences between the PDFs (see figure 
The estimates for the three different cutoff values are consistent within the errors, indicating a weak dependence on f c . 



3.3 Estimated parameters from Q1422 

In the top panel of figure H we plot the normalized flux of the two spectra of Q1422 obtained independently by Rauch et 
al. (solid line) and Songaila & Cowie (dotted). There is a good general agreement between the two. As we will see later, 
the discrepancies caused by the different continuum fits have a moderate effect on the estimated parameters. In the bottom 
panel of figure @ we plot the recovered local optical depth in redshift space. The results for both spectra agree very well for 
t Si 3. Both curves have a cutoff in saturated regions where the flux goes to zero. The differences at f > 3 will not affect our 
estimates of B and C because of the truncation at smaller values of f . 

Figure (0) shows fP[ln(f )] for the recovered f from the Rauch et. al. spectrum, in real (solid line) and in redshift space 
(dashed). The dotted line is for the log- normal with parameters estimated from f in real space. There is reasonable agreement 
between the recovered PDF and the log-normal for low optical depth implying that the true PDF cannot be very different 
from a log-normal for low densities. 

In table (|3j) we list the estimated values of C and B and the value of A inferred for a — 1.7. Because of the strong redshift 
dependence of A (eq. |4J) , we present results from the first half, the second half and the total length of both spectra. We have 
used a variable Doppler parameter b = 60 f 1-c * ' 2 with a — 1.7. A value of 60 = 15km/s matches the \ 2 criterion described in 
the previous section. As discussed in section 3.2. our inversion procedure works well for cut-off values r c up to 2. For t c = 2 we 
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Figure 5. The normalized fluxes (top) and recovered local optical depths in rcdshift space (bottom) for a fraction of the two spectra 
of Q1422. The solid and dotted lines correspond to the spectra observed and continuum fitted by Rauch et. al. and Songaila & Cowie, 
respectively. 
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Figure 6. Curves of f P[ln(f )] for the local optical depth in real (solid) and redshift space (dashed) recovered from the Rauch et. al. 
spectrum. The dotted curve corresponds to the log-normal PDF obtained from the estimated parameters (C, B) = (1.8, —1.61). The error 
bars are taken from the figure 2. 



Table 3. Estimated parameters from the Rauch et al (R) and Songaila & Cowie (SC) observations of the spectrum of ql422+231. Listed 
are (C B, A) where A are computed from C and B using ln(A) = B + C 2 /2/a and assuming a = 1.7. The uncertainty in these parameters 
is estimated as the maximum total error from the 4 sets of the 70 mock spectra (see text). The absolute lcr on Cand B are 0.27 and 0.4. 







fc = 1.0 


r c = 1.5 


f c = 2.0 






(C,B,A) 


(C,B,A) 


(C,B,A) 


Low redshift half 










realspace, R 




(1.83,-1.59,0.54) 


(1.76,-1.74,0.43) 


(1.70,-1.82,0.38) 


real space, SC 




(2.21,-1.09,1.44) 


(2.09,-1.78,0.60) 


(2.05,-1.87,0.53) 


redshift space, 


R 


(1.82,-1.91,0.39) 


(1.79,-1.90,0.38) 


(1.71,-1.97,0.32) 


redshift space, 


SC 


(2.35,-1.71,0.92) 


(2.22,-2.03,0.56) 


(2. 15,-2.132, .46) 


High redshift half 










real space, R 




(1.89,-0.89,1.17) 


(1.89,-1.32,0.76) 


(1.87,-1.35,0.72) 


real space, SC 




(2.06,-0.63,1.86) 


(2.12,-1.26,1.06) 


(2.10,-1.34,0.95) 


redshift space, 


R 


(1.99,-1.35,0.83) 


(1.97,-1.44,0.74) 


(1.89,-1.62,0.56) 


redshift space, 


SC 


(2.29,-0.95,1.81) 


(2.26,-1.55,0.95) 


(2.126,-1.84,0.60) 


All 

real space, R 




(1.89,-1.26,0.81) 


(1.84,-1.55,0.57) 


(1.80,-1.61,0.51) 


real space, SC 




(2.14,-0.84,1.66) 


(2.13,-1.52,0.83) 


(2.09,-1.61,0.72) 


redshift space, 


R 


(1.92,-1.66,0.56) 


(1.89,-1.70,0.52) 


(1.81,-1.82,0.42) 


redshift space, 


SC 


(2.35,-1.33,1.34) 


(2.26,-1.79,0.75) 


(2.14,-1.98,0.53) 



obtain (C,B) = (1.8, —1.61) and (C,B) = (2.09, —1.61) for the Rauch et al and the Songaila & Cowie spectrum, respectively. 
As discussed in section 3.2. the value of B is expected to be biased low. 

To assign errors to the estimates of C and B we draw on the analysis of the mock spectra in the previous section. We 
take the maximum error of all sets of mock spectra calculated for half the length of the observed spectrum as a conservative 
error estimate. This gives ac = 0.27. and erg = 0.4, respectively. The values of C and B derived from the upper and lower 
half of the spectrum are consistent within our estimated errors if we account for the expected redshift evolution of B. The 
redshift evolution of the clustering amplitude and a and therefore that of C is expected to be small. The difference between 
the estimates obtained for the Rauch et al and the Songaila & Cowie spectrum are also consistent within the estimated errors. 

For a given value of a we can now infer values of a and A from our estimates of C and B. The values of A given in table 
(H) are derived using 

A = exp[B + C 2 /2a\. (12) 

Since B is systematically biased low, the values of A in the table should be considered as lower limits. Assuming that the true A 
is biased low by a factor of 2 as suggested by our analysis of the mock spectra with /i3 = 0.4 we estimate A(z w 3.25) « 1.2±0.7 
(mnA = >/(<x| + acC/a) « 0.65 for a = 1.7). 

We now discuss the inferred value of a =< 8 2 > 1//2 . For a PDF described by an Edgeworth expansion (equation hll), a 
and /i2 =< (log[l + <5]) 2 > ' 2 are related by 

a 2 - (l + ^2^3/3!) exp (/4) - 1. (13) 

We have demonstrated in the previous section section that the estimates of C are unbiased even for true PDFs with rather 
large skewness parameter ^3. The value of (12 — C/a is therefore also unbiased. Hence by assuming /13 = we somewhat 
overestimate a. For the values of /12 and (13 measured in numerical simulations (NH99) we expect this bias to be less than 
15%, significantly smaller than the uncertainty in the value of a. We therefore use the relation 

a 2 = exp [(C/a) 2 } ~ 1. (14) 

Note that the small bias introduced here will also somewhat lower our estimates for fi in the next section. Possible values of 
a range between 1.56 and 2 but a is probably about 1.7 at z ~ 3 for plausible reionization histories (Hui & Gnedin 1997, 
Abel & Haehnelt 1999). The estimated value of C = 1.8 ± 0.27 for the Rauch et al. spectrum translates into a(a — 2.0) = 
l.ll±0.2,a(a = 1.7) = 1.44±0.3, and a(a — 1.56) = 1.67±0.4. The values for the Songaila & Cowie spectrum are 20 percent 
higher. 
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Figure 7. The non-linear rms value of density fluctuations as a function of the effective Jeans length scale. The solid, dotted and dashed 
curves correspond to T = 0.21, 0.1 and 0.5, respectively. The three points with error bars are the values estimated from the Rauch ct al 
spectrum of Q1422 for a = 1.5, 1.7 and 2.0 (from high to low). The abscissa of the data points are arbitrary and the errors are la and 
include the error due to cosmic variance and the inversion. 



4 COMPARISON WITH THE MASS FLUCTUATIONS AMPLITUDE PREDICTED BY CDM-LIKE 
COSMOGONIES: IMPLICATIONS FOR fi 

There is a whole family of CDM-like power spectra which is consistent with a variety of observational constraints on the DM 
clustering amplitude. If we normalize these power spectra to match the local space density of galaxy clusters (e.g., Eke, Cole 
& Frenk 1996) then, for a given fl, we can predict the fluctuation amplitude on the scale and redshift relevant to the estimate 
of a in the previous section. In the following we obtain constraints on Q by comparing our measurement of a with such model 
predictions. 

We start with the linear mass power spectrum as given by Bond, Eftsathiou, Frenk & White (1993) 

Ak 



P(k) = 



(15) 



[1 + (ak + (6fc) 3 /2 + (cfc) 2 )"] 27 " ' 

where the wavenumber k is in units of l/(/i Mpc), A is a normalizing factor, a = 6.4r~ 1 /i~ Mpc, b = 3.0T ft Mpc, 
c = 1.7r _1 /i _1 Mpc, and v — 1.13. For standard cold dark matter models the constant F is identified with fih. The space 
density of present-day galaxy clusters is reproduced if A is chosen such that the linear rms value, as, of the density fluctuations 
in a sphere of radius 8/i _1 ft _1 Mpc is given by 

-0.52+0.130 _ _ tn KO j_ n nA ^ ,-,-0.46+0. m 



as = (0.52 ± 0.04) fi" 



(78 = (0.52 ± 0.04) fr 



(16) 



for flat and open models, respectively (Eke, Cole & Frenk 1996). We then use an analytic transformation from the linear to 
the non-linear power spectrum (Peacock & Dodds 1996; see also Jain, Mo & White 1995). 

In section 2 we introduced an effective Jeans scale x^ such that the fluctuation amplitude of the DM density field smoothed 
on this scale can be computed from the non-linear power spectrum as 

-2 



c^O.,) = 



(27T) 



dfc 3 P nl (fc) 



1 + 



2tt 



(17) 



Even though x } should be related to the linear Jeans scale, the exact choice of x } is not obvious and we have left x } 
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Figure 8. Contours of x = 1 (solid) and x = 3 (dotted) as functions of f2 and x 3 for two values of a for open and flat models, 
as indicated in each panel. The x 2 function quantifies the deviation of the measured clustering amplitude from that predicted by the 
r = 0.21 model at 2 = 3.25. The error includes cosmic variance, the error due to the inversion and the uncertainty in the normalization 
of the linear power spectrum. 

as a free parameter. The nonlinear a depends then on the shape parameter F, the background density, fl, the cosmological 
constant, A, and the effective Jeans scale, a; T . In figure (hj) we show the model prediction for a at z = 3.25 as a function of the 
effective Jeans scale Xj for different cosmological parameters. The value of a inferred from the Rauch et al. spectrum is shown 
as the open symbols with error bars for a = 1.5, 1.7 and 2.0. Three factors contribute to the rather strong Q, dependence of 
the predicted a: the power spectrum normalization, the time evolution of the fluctuation amplitude and the relation between 
length and velocity scale. There is a similar but somewhat weaker dependence on A. 

For a given a, we define the function \ 2 — [C/a — ln(l + a 2 )/2\ /cr 2 rr where the <r crr is the total la error including the 
uncertainty in the power spectrum normalization Qlq). Figure (H) is a contour plot of x as function of Q and a;, for flat 
(SI + A = 1) and open models (A = 0). If we assume V — 0.21 as inferred from the shape of the power spectrum of the observed 
galaxy distribution at z = (Peacock & Dodds 1994) our estimate for S7 can be expressed in a convenient parametric form 
as follows: S7 = 0.61 (a/1.7) 1 ' 3 (ij/0.62)" ' 6 for a flat universe, and Q = 0.91 (q/1.7) 1 ' 3 (x J /0.62)~ ' 7 for an open universe 
without a cosmological constant. The errors on Q are 25%. The inferred values from the Songaila & Cowie spectrum have the 
same parametric form but are lower by 15%. 



5 DISCUSSION 

5.1 Comparison to previous work 



The value of A. 

Rauch et al. have used several HIRES spectra including the one of Q1422+231 used here to infer the value of fi — 
(flbarft 2 /0.0125) 2 (ii"(2)/100kms~ 1 Mpc _1 )" 1 (r phot /10~ 12 s _1 )~ 1 by comparing the flux decrement distribution of the observed 
spectra to those of mock spectra calculated from hydrodynamical simulations. This analysis was basically equivalent to a 
determination of the optical depth constant A. The two quantities are related as 



A(z) « 0.71 n{z) (f/1.5 x 10 4 K)-°- 7 ((1 + z)/4) 6 



(18) 
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Rauch et al obtained fj,{z ~ 3.25) ss 0.6 and 1.5 for mock spectra produced from a SCDM model (Hernquist et al. 1996) and 
ACDM model (Cen et al. 1994) . The temperatures at z = 2 were t = 5.6 x 10 3 K and f = l.lx 10 4 K. This corresponds to 
A(z — 3.25) ~ 1.2 and and 2.0, respectively. No correction for the temperature evolution in the simulations has been made. 
This is in reasonable agreement with our estimate of A ~ 1.2 ±0.7 but we would like to stress that our value includes a rather 
large correction of a factor two for the bias in B which we have inferred from our analysis of mock spectra produced with a 
PDF with p3 — 0.4. We therefore do not think that we can currently improve the constraints of Rauch et al. on the baryon 
density and the flux of ionizing photons. The fact that we obtain a consistent result with a completely different procedure 
is nevertheless very encouraging. Our value is larger than the A(z ~ 3.25) ~ 0.4 in the fiducial model of Bi & Davidsen 
(a = 1.77). 

The value of a 

Gnedin and Hui (1998) used mock spectra calculated from DM simulations to investigate the influence of the fluctuation 
amplitude on the slope of the column density distribution of QSO absorption lines. They got reasonable agreement with the 
observed slope if the linear fluctuation amplitude at half the linear Jeans scale for T = 10 4 K was about 2. This is larger 
than our value for the non-linear a. This indicates that the effective Jeans scale is larger than half the linear Jeans scale for 
T = 10 4 K. Bi& Davidsen have assumed a = 1.95 in their fiducial model significantly larger than our value a(a = 1.77) = 1.34 
for the same value of a. 

Croft et al. (1998, 1999) have compared the flux power spectrum of observed absorption spectra to infer the DM power 
spectrum on scales between 2/i -1 Mpc and 12/i _1 Mpc (comoving, fi = 1 ) by comparison with mock spectra of DM simulations 
which were calibrated using hydrodynamical simulations. Note that Croft et al. define their length scales to be a factor 2-7T 
larger than we do in our definition of the Jeans length. Croft et al obtain A = k 3 P(k)/2ir 2 ~ 0.57 ± 0.3 at z ~ 2.5, where 
A is the rms fluctuation amplitude of the DM density contributed by a logarithmic interval of k. This measurement was at 
a scale which corresponds to x « 12(f2/i 2 /0.125)~°' 5 Mpc. Weinberg et al. (1999) then used the constraints on the DM power 
spectrum obtained by Croft et al. (1999) to constrain the value of fl with similar assumptions for the normalization and the 
shape of the power spectrum as we have made here. For T — 0.2 Weinberg et al. 1999 obtain Q ~ 0.34 ±0.1 (A = 1 — fl) and 
Q ~ 0.46 ±0.1 (A = 0). We obtain somewhat higher values unless the effective Jeans scale is larger than the linear Jeans scale 
(for a temperature of 15000 K) by about a factor two. 

5.2 Future work 

There are a number of possible improvements on the analysis presented here. As discussed in section 3.4 the true PDF of the 
density is expected to be somewhat positively skewed relative to log-normal, but we were not able to constrain the skewness 
parameter, /13 . To make progress on this issue an improved model of the PDF which e.g. relates j-13 to other parameters would 
be needed. There are attempts to determine the value of a from the correlation of column density and the lower cut-off of the 
Doppler parameter distribution (Schaye et al. 1999). Such a determination will obviously improve our constraints on A and 
a. Currently we probe the density PDF up to a rather low overdensity of S ~ 2. This could be extended to higher densities 
by incorporating higher order Lyman series line into the inversion scheme. This will tighten the constraints on the density 
PDF and also extend the redshift range of each spectrum available for the inversion. The main source of uncertainty in our 
constraints on O is the effective Jeans lenght scale x T . A calibration of x , by numerical simulations should help to tighten the 
constraints on Q. It should also be possible to extract information on the shape of the power spectrum on small scales from 
the correlation function of the truncated real-space optical depth. We estimated the error of a due to cosmic variance to be 
about 20% for one spectrum like Q1422. This error could be reduced to 5% by applying the method to 20 quasar spectra. 
Note that the error estimate was based on mock spectra generated with a SCDM power spectrum. It might be different for 
other power spectra. It would obviously also be interesting to probe the redshift evolution of the clustering amplitude. 

6 CONCLUSIONS 

We have recovered the real space optical depth distribution along the LOS to the QSO Q1422±231 from two different HIRES 
spectra of Q1422+231 using the inversion technique proposed by Nusser & Haehnelt (1999). Despite the two spectra having 
different continuum fits, the difference in the recovered optical depth distribution is small. The recovered optical depth 
distribution is similar to that expected for a log-normal PDF of the density. 

We have used a new technique to estimate the fluctuation amplitude of the gas density and the optical depth constant 
from the truncated optical depth distribution. Assuming that the PDF of the gas density is log-normal we showed that the 
moments of the distribution depend on two parameters B — In [.4/(1 + c) 2 ] and C — ayln [a 2 ± 1]. This demonstrates that 
there are fundamental degeneracies between fluctuation amplitude of the gas, the power law index a, and the optical depth 
constant A for a log-normal PDF. These degeneracies will be only very weakly broken if the true PDF is somewhat positively 
skewed relative to log-normal as expected from numerical simulations. 
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Using mock spectra we have shown that robust estimates for the parameters B and C can be obtained. The mean of the 
recovered C is unbiased with respect to the true value. The value of B however, is systematically biased low if the true PDF 
is positively skewed relative to log-normal. Our estimates of B and C translate into A ~ 1.2 ± 0.7 and a ~ 1.44 ± 0.27 if 
we assume a = 1.7, the most likely value for plausible reionization histories. The error estimates come from our analysis of 
the mock spectra. If the IGM indeed traces the DM in the way suggested by simulations this is the first measurement of the 
fluctuation amplitude on an effective Jeans scale which corresponds to about two hundred kpc(comoving). The fluctuation 
amplitude of the DM density on such small scales is of particular interest for the question when the epoch of reionization 
occurs. Our rather low estimate for the value of a strongly argues for a late epoch of reionization. 

Finally, we have compared our measurement of the fluctuation amplitude to predictions by CDM-like cosmogonies with 
a power spectrum shape parameter suggested by the clustering behavior of galaxies, and normalized to match the observed 
space density of present-day clusters of galaxies. The fluctuation amplitude of such models depends rather strongly on Omega. 
We obtain fi = 0.61 (a/1.7) 1 - 3 (z T /0.62)~ ' 6 (A = 1 - fi) and fi = 0.91 (a/1.7) 1 - 3 (aij/0.62) -0 - 7 (A = 0). If we assume that 
the effective Jeans scale is close to the linear Jeans scale (at 15000K) then the measured fluctuation amplitude is larger than 
that predicted for an (Q = 0.3, A = 0.7) Universe by a factor of about two. This argues either for a larger value of Q and/or 
a large value of the effective Jeans length x } . 
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